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With the general theory of standard quantum process tomography, we shall develop a scheme to 
decide an arbitrary matrix element of \, which is in the Choi matrix representation, in a scalable 
way. 
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The characterization of the evolution of a quantum sys- 
tem is one of the main tasks to accomplish to achieve 
quantum information processing. A general class of 
methods, which have been developed in quantum infor- 
mation theory to accomplish this task is known as quan- 
tum process tomography (QPT)- for a review of quan- 
tum tomography, see Refs. [1 — 3]. Under general as- 
sumptions, the evolution of quantum system can be rep- 
resented by a linear, completely-positive map which can 
be written as 



(1) 



with £) m (J5 m )tJ5 m = I. Using a fixed set of operators 
{Ei}, which form a basis for the set of operators on the 
state space, so that E m = £\ e™Ei, one may introduce 
the x matrix representation with Xij = Em e T( e T)* f° r 
its matrix elements, and rewrite (1) in the way like 



e(p)=J2x^E l (p)El 



(2) 



The task of QPT can be viewed as to decide \ y i a 
different protocols. Standard quantum process tomog- 
raphy (SQPT) was the first method proposed [1,4,5]. 
It involves preparing a set of input {pj} and measur- 
ing its output s(pj) via the quantum state tomogra- 
phy. With the denotations that s(pj) = J2k^3 k P k ano ^ 
EmPjEjn — J2j Pij n Pj> where /3™ n are complex numbers 
which can be determined by standard algorithms from 
linear algebra given the {E m } operators and {pj} oper- 
ators [1], one may get the relation 
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(3) 



where Xmn can be determined by the given parameters 
/3£? n and the data of known from quantum state to- 
mography. Besides SQPT, other methods like the ancilla- 
assisted quantum process tomography have also been 
proposed [6]. 

Usually, the complete characterization of \ matrix is 
a non-scalable task: For the N d-levels system, there are 
about d AN elements to be decided. Recently, a series of 
works have demonstrated that it is possible to extract 
partial but nevertheless relevant information about the 
quantum process in an efficient and scalable way [7-13]. 



These approaches share an essential feature that: They 
are based on the idea that the tomography of a quan- 
tum map can be efficiently performed by studying cer- 
tain properties of a twirling of such a map. Another 
method, the so-called direct characterization of quantum 
dynamics (DCQD), was also constructed for use in par- 
tial characterization of quantum dynamics [14-16]. 

Here, we shall focus on the topic whether the general 
theorem of SQPT developed in [1] can be applied for par- 
tial QPT. The case, where the \ is defined by the Choi 
operators, should be discussed in present work. With the 
general theory of SQPT in [1], we shall develop a scheme 
to decide an arbitrary \ matrix element, which certainly 
carries some partial information about the quantum pro- 
cess [13], in a scalable way. We shall show that: For the 
X in the Choi matrix representation, to decide its diago- 
nal element one need just a single measurement while to 
decide its off-diagonal we require about sixteen measure- 
ments. This result is independent of the actual dimension 
of the system. Our scheme also represents a new kind of 
method of partial QPT without of ancilla. In following 
argument, we shall to show how to construct our scheme 
in detail. 

The general theorem.- With {|a)} a =o,i d-i the ba- 
sis of the D-dimensional system, we introduce the Choi 
matrix 



E ab = \a){b\. 



(4) 



This set of operators, {-E a b}a. h=o,i D-i, forms an or- 
thogonal basis with the algebra that Tr(£' ( luEd) = 
S ac Sbd- Any D x D matrix, say, M, can be expanded 
with it, M — y^~Q MatjEab, where the expanding co- 
efficients M a b = Tr{E\ b M) are just the matrix ele- 
ments of M. Let each Kraus-operator in {E m } to be ex- 
panded with the conjugated Choi matrices in (6), E rn — 
EeJ= e?f E ef with el) = Tr[El f E m ] and (E m ^ = 

Ylg h=o( e g%)* (EghY, we can define the x m the Choi ma- 
trix representation with its elements to be 



Xef;gh = V < 



efV^gh 



(5) 



We call Xef;gh the diagonal matrix element if e = g and 
/ = h, else, we call it the off-diagonal matrix element. 

The relation in (3), which is one of the main results 
of the general theorem of SQPT [1], puts no constraints 
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on our choices to define A and x matrices. In present 
work, we use {E ab } to describe the input and e(E ab ) 
for its output. Let e(E ab ) be expanded with the set of 
Choi operators in (4), e(E ab ) = Y^c,d=a ^ab;cdE cd with its 
elements to be 

X ab:cd = Tr[El d e{E ab )]. (6) 

With the above definitions in hands, it can be shown that 
there exists a one-to-one mapping between the A matrix 
element and the \ matrix element, 

X a b;cd = Xca;db- (7) 

This result can be proved by the general theory of SQPT. 
From the expressions of X ab[C d and Xef;gh m above, there 
should be 

D-l 

X a b;cd = Pll-cdXef-.gh (8) 

e,f,g,h=0 

where the transformation matrix (3 has its elements to be 

P:t:$=TT(El d E ef E ab El h ). (9) 

Each matrix element of [3 can be directly calculated, 
& = Tr[M)( C |) • |e)(/| • \a)(b\ ■ (\h)(g\)} = 5 ec 8 fa 8 gd S hb . 
With columns indexed by ef; gh and rows by ab; cd, (3 
can be expressed as a D A x D A matrix. In each row and 
each column of it, there is just one non-zero matrix ele- 
ments (with the value of 1). Furthermore, one may verify 
that the determinant of (3 equals 1 while its inverse is the 
transpose of itself, 

det(/3) = 1,/3 _1 = (10) 

According to the famous Cramer's rule [17]: If Ax = b 
is a system of n linear equations in n unknowns such 
that det(^4) ^ 0, then the system has a unique solutions. 
Equation (8) can be viewed as such a set of linear equa- 
tions with: A -» /3, x -> x, an d b -> A. From (10), 
we conclude that x IS unambiguously determined by A, 
X = /3 T X, since det(/3) = 1. Finally, using the known re- 
sult of Pal'cd f° r (8) j we shall get the one-to-one mapping 
described by (7). 

As a direct application of the one-to-one mapping 
in (8), we find that the diagonal matrix element 
Xab;ab(^bb;aa) is physical meaningful: it represents the 
probability of the transition from the initial a-th level to 
the final b-th level, 

Xab;ab = (b\e(\a)(a\)\b). (11) 

In other words, the classical measurement of the transi- 
tion probability of the D-levels atom can be also viewed 
as to decide the diagonal x matrix element in the scheme 
of SQPT. With a simple reasoning, we have 

D-l D-l 

Trx = £ Xab-,ab = ]T Tr[e(|a) (a|)] . (12) 

a, 6=0 a=0 



For the trace-preserving cases, Trx = D. 

Choi matrix SQPT.-As we have shown, the x ma- 
trix in (5) should be given if the A matrix in (6) has 
been decided. However, the A matrix elements, X ab -.cd — 
Tr[E\ d e{E ab )], can not be directly measured if one of the 
operators, E\ d and E ab , is non-Hermitian. This prob- 
lem can be solved with the following protocol: At first, 
one may introduce a set of linearly independent states, 
{|5' m )} m=1 D 2, and expand each E ab as 

D 2 

E ab =Y,rt\*m){*m\ (13) 

m—1 

where the coefficients are known. Now, e(E ab ) — 

D 2 

J2 m =i r m s (\^m)(^m\) since that the operation of e is 
linear. Then, giving a set of linear independent Hermi- 

tian operators {O n }n=i,...,D 2 (O n = O^), we rewrite & cd 
as 

D 2 

Eld = J2 s n 6 - (14) 

n=l 

with known parameters s^f. Because that perform- 
ing trace is also a linear operation, there should be 
TT[E cd (...)} = E^i4 d Tr[0„(...)]. Finally, the way of 
measuring X ab . cd (xca-db) is clear, 

D 2 

\ab;cd= £ rt S c n d Tr[6 n e(\^ m )^ m \)]. (15) 

m,n—l 

The configuration space of the measurements is known, 

M : {Tr[d„£(|V> m )(^ ro |)]} TO ,„=i,2,...,.D2- 

To perform the complete SQPT, we shall carry out all 
D 2 x D 2 measurements in M and decide each A(x) ma- 
trix element one by one according to their known co- 
efficients, rf b and Sj d , here. However, if only a selected 
X a b;cd(Xca;db) is to be measured, one may just perform the 
measurements satisfying the constraint that r^ b ■ s^ ^ 0. 

In fact, it has been shown that each E ab can always be 
expanded with four pure states [1]. Introducing 

\ab, +) = ^(\a) +\b)), \ab, -) = ^(\a) + i\b)), (16) 

for a < b, each non-Hermitian E ab can be expanded with 
Eab = \ab,+){ab,+\ +i\ab, -){ab, -| 

-^(\a)(a\ + \b)(b\). (17) 

For the case a > b, E ab can be directly de- 
rived with E ab = El a . Formally, for a ^ b, 
we have E ab = £? =1 rf^f ){^f | with V|Vf) € 
{|a), \b), \ab, ±)}. In the similar way, we expand Er cd 
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with E 



Id = Ej=i sf\^f)(ipf\ by requiring V|^ cd ) G 
{|c), |cd, ±)}. Now, equation (15) is simplified into 



A a6;cd - J2 rfsf{ipf\e(\ipfm b \Mf). (18) 

ij = l 

To decide the off-diagonal matrix element of x, we need 
sixteen measurements by taking IV'f 6 ) for the input state 
and measuring its output state e(|V , f b )(' i /'f f 'l) with the 
projective operator \ipj)(ipj\- 

It can be seen that the number of all the possible input 
states, \ipf) with 0<a<6<D-land j=l, 2, 3, 4, is 
limited to be D 2 : The basis vector \a) has D terms, while 
-^2(|a) + \b)) and ^(|a) + i\b)) has the same number of 
\D(D — 1), respectively For convenience, we use S to 
denote all these states, 



S:{\a),^(\a) 



\b)),%\a) 



i\b))} 



0<a<b<D-l- 



For the trace-preserving cases, only D 2 (D 2 — 1) measure- 
ments are independent. An interpretation for it is like 
this: the following D measurements, {a\e{\tp°; b ){'4>'} b \)\a) 
with a=0, 1,..., D-l, should be performed for a given in- 
put state \tpi b }- With J2a=o l a )H = we can a l wavs 
leave (D— l\e(\ijj° ,b } {tpf b \)\D — 1) unmeasured and decide 
its value according to 



(D-l\e(Wf)(i>f\)\D-l) = l 



D-2 

E 

a=0 



(a\s(\ipf)(W»\)\a). 



This fact is consentient with our general analysis about 
the x matrix: For the D-dimensional system, it has 
D 2 (D 2 — 1) independent real parameters and D 2 addi- 
tional constraints for trace. 

The N d-levels system.- The way of performing partial 
QPT is important for the the N d-levels system specially 
when d is a large number. For such cases, to performing 
the complete QPT becomes a non-scalable task since the 
number of the required experiments is exponentially in- 
creased with N. Our general protocol also holds for the N 
d-levels system by taking it as a special case of D = d N . 
With 



«i 



|aj) ® • • • ® \a N ) 



(19) 



where {|aj)} aj= o.i,...d-i is the basis of the j-th d- 
dimensional subsystem, we could define a relation be- 
tween the single index with its corresponding string of 
local indices, say, a — > aia 2 ■ ■ ■ ajv, b — > &i& 2 ■ ■ ■ &iv, etc.. 
If the set of product basis {EL/Li ®\ a j)} 1S to be used in- 
stead of {\a)} used above, all the equations in the product 
basis can be easily given by substituting each single index 
with its corresponding string of local indices. For exam- 
ple, the one-to-one mapping in (7) with the local indices 
should be, 



For the two qubits case, with |+) = ^(|0) + |1)) and 

|— ) = ^(|0) +i\l)), S contains following product states, 
|0)|0),|0)|1),|1)|0),|1)|1),|0)|±),|1)|±),|±)|0),|±)|1), 
and the four Bell-type states, ^(|00) + |11)), 

^(|00) + i|ll)), ^(|01) + |10», and ^(|01) + i|10». 
For the more general N d-levels case with the product 
basis, S contains a series of product states and M-parties 
maximally entangled states with M arranged from 2 to 
N. Suppose we have the freedoms of choosing one of the 
followings two mappings, |aj) — > | ti)> \h) — > I 4-*) or 
\o-i) — > I ti)j ~^ I 4-i); f°r each site, by renumbering 
the sequence of all the N sites, we can write \ab,±) in 
(16) with 



N 



N 



\ab,+) = \G M ) [] ®|t,),K-)HG M 



j=M+l 



II ®lt;> 

j=M+l 



where the two Greenberger-Horne-Zeilinger type states, 
\G M ) = \)® M + | |)® M ) and \G M ) = ^(| f)® M + 
i\ |)® M ), are the maximally entangled states among M- 
parties. 

Discussion- With the general theory of SQPT in [1], 
we have developed a scheme of performing partial SQPT 
for x i n the Choi matrix representation. Besides the fact 
that the set of Choi operators is a convenient basis for an 
arbitrary D-dimensional system, x matrix in this repre- 
sentation is shown to be physical meaningful in the sense 
that its diagonal matrix element is just the transition 
probability of the system. For the N qubits case, it can 
be shown that x° (m the Choi matrix representation ) 
and x P (m the Pauli matrix representation), are equiva- 
lent with each other. An interpretation is like this: Let's 
at first consider the single qubit case, it can be seen that 
the set of Pauli operators is related with the Choi oper- 
ators by a unitary transformation U, 
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where we define E a b = \/2E a b. One may easily veri- 
fied that l/W = I4 and get the relation x P = Ux C U^ . 
For the N qubits case, if x P is expanded with the set 
of Pauli operators {IlfcLi ®v k } while x° i s defined with 
{JlfeLi ®E k }, the relation X P = Ux°U^ still holds with 
U = rifcLi U k ■ From above discussion, it's hard for us to 
claim that one representation has privilege over the other 
representation for defining x matrix. 

From equation (15), it should be noted that the ways 
of performing complete SQPT are not limited, any set of 
linearly independent states {l^m)} with another set of 
linearly independent Hermitian operators {O n } can be 
applied for this task. Among all the possible ways, our 
protocol in (18) works for deciding an arbitrary x ma_ 
trix element with a fixed number of measurements. If 
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the complete SQPT should be performed, one may use 
the following scheme: Let {\<pj}j—i,...,d 2 the set of lin- 
early independent states for the the k — th subsystem, 
the input states \^ m ) in (13) may be defined as a pure 
product states, say, l^m) = Ylk^i^l^)' Considering 
that f§ = I d and {t k - }j=i,...,d2_i, where is the gen- 
erator of SU(d), forms a orthogonal basis of the k-th 
subspace, the operators O n in (14) can be chosen to be 
O n = IlfcLi ®r fe , the complete SQPT performed in this 
way has the unique property that it does not require any 
two-bodies or multi-bodies interaction. 

There are two related problems still unsolved with the 
present work. At first, if all the experiment data are used 
to reconstruct the whole \ matrix, how to keep the pos- 
itivity of \ i s n °t given since that the partial SQPT is 
concerned here. Second, our protocol of partial SQPT is 



limited to the case where \ is m the Choi matrix repre- 
sentation, it is still an open problem whether the general 
theory of SQPT in [1] can be also applied for the partial 
QPT when \ is m other physical representations. 

Finally, let's make a short summary for our work. With 
the general theorem of SQPT, we developed a method to 
estimate an arbitrary \ matrix element in the Choi ma- 
trix representation. Our observation is that: To decide 
the diagonal matrix element, one just needs a single mea- 
surement; To decide an arbitrary off-diagonal matrix el- 
ement, we should carry out sixteen measurements. This 
result is independent of the actual dimension of the sys- 
tem. Compared with the known methods of partial QPT, 
our scheme does not require any additional resources. It 
can be applied for the case where a clean ancilla system 
is not available. 
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